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Abstract 

The recent experimental indications of density dependence in the pion decay con- 
stant /* and u meson mass m* and the discovery of S^SHS) and other strange 
"nuggets" are providing a strong support for kaon condensation in dense hadronic 
matter, thereby re-kindling the interest in the issue of the critical stable mass of neu- 
tron stars. The density-dependent quantities provide increases in the vector mean 
fields mediated by p and w-meson exchange which increase by a factor ~ 1.56 the 
Weinberg- Tomozawa term in kaon-nucleon interactions which accounts for ~ half 
of the binding energy of the K~ meson in dense matter. Furthermore lattice gauge 
calculations have pinned down the value of KN sigma term, T^kn, the explicit chi- 
ral symmetry breaking in the strangeness sector. The partial rotation out of this 
explicit breaking provides the other ~ half of the K~ binding energy. The net re- 
sult is to confirm the work of Thorsson et al. that strangeness condensation takes 
place at u = u/uq ~ 3, where no is nuclear matter density, in neutron stars. We 
suggest in this article that a support for this scenario is provided by the recent 
experiments of Suzuki et al. who found tightly bound strangeness nuggets. The 
strangeness nuggets discovered in the experiments involve approximately the same 
ratio of nucleons to K~ meson as in the center region of neutron stars. But whereas 
the latter can be described by mean fields, in which the medium effects are substan- 
tially more attractive than in the finite system, the binding in neutron star matter 
should be substantially (say, ~ 20%) greater than that in the strangeness nugget. 
This would strengthen the argument by Brown and Bethe that the accompanying 
softening in the equation of state should limit the maximum neutron star mass to 
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M^g x ~ 1.5M Q . This low M^ x ~ 1.5M© has major consequences in astrophysics, 
especially for the merging rate of compact stellar objects. 



1 Introduction 

Based on the detailed numerical calculations of Thorsson et al. [1], Brown and 
Bethe [2] claimed that the result of kaon condensation setting in at a density 
n ~ 3n , where n is nuclear matter density, gives a sufficient softening of 
the equation of state in neutron star matter to limit the maximum possible 
neutron star mass to ~ 1.5M . In the analysis of Supernova 1987A, which 
Bethe and Brown considered to have evolved into a black hole, these authors 
[3] were able to establish a maximum mass for the compact object of 1.56M , 
based on the observed amount of nickel produced in the supernova explosion. 
The M™g x ~ 1.5M was consistent with 1987A evolving into a low- mass black 
hole. Furthermore, the ZAMS (Zero Age Main Sequence) 18M progenitor of 
1987A was calculated to have an Fe core mass of ~ 1.5M . Brown et al.[4,5] 
argued that this Fe core mass was about the same as the compact object mass 
on the basis of calculations by Woosley showing that fallback in the supernova 
explosion compensated for the increased binding energy of the compact object. 

We claim that the above astrophysical phenomenon involving neutron star 
masses is intricately connected with the softness in the equation of state com- 
ing when the electrons in the neutron star can shed their high degeneracy 
energy by changing into i^~-mesons, in a zero momentum Bose condensate. 
This is a particle and nuclear physics problem which in our view must incor- 
porate the vacuum structure of dense matter involving chiral symmetry and 
hence in-medium properties of hadronic masses. Our focus here will therefore 
be on particle/nuclear physics aspect. We shall make a brief summary of the 
astrophysical consideration in a section below but leave the detailed discussion 
to a longer, technical paper [6]. 

The standard chiral perturbation approaches predominantly employed in the 
field are anchored on chiral Lagrangians whose parameters are determined in 
the matter-free vacuum (that we shall refer to as "zero- vacuum" ) and compu- 
tations are done by perturbation around the zero-vacuum in terms of a set of 
presumed small expansion parameters. Our approach will differ from them in 
a crucial way. Ours will be a mean-field approach with a Lagrangian defined in 
a sliding vacuum [7]. The basic idea is similar to Landau Fermi-liquid theory 
for many-body systems as applied to nuclear matter [8]. Here one performs 
"double decimation" [7]. The crucial point is that the mean field approxima- 
tion with an effective chiral Lagrangian with scaling parameters is equivalent 
to Fermi-liquid fixed point theory. It is not known how to obtain the effective 
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Lagrangian with "sliding vacua" from an effective Lagrangian defined in the 
zero- vacuum just as it is unknown how to derive from a fundamental chiral La- 
grangian the four- Fermi interaction Lagranian effective near the Fermi surface 
that via Wilsonian RGE gives rise to Fermi-liquid fixed point theory [9]. In the 
same vein, our starting Lagrangian that is defined with "sliding vacua" is not 
derived in any systematic way from a zero- vacuum chiral Lagrangian. It will be 
the same Lagrangian that gives rise to the Fermi-liquid fixed point theory, here 
extended to three flavors. We will employ mean field approximation with this 
Lagrangian. Kaon condensation will emerge as instability against strangeness 
condensation like superconductivity resulting from an unstable Fermi surface 
triggered by the Cooper pairing. Since in this way of approaching kaon con- 
densation one is approaching from below the QCD phase transition, namely, 
chiral restoration, we will call this "bottom-up" approach. 

Since kaon condensation must take place before, and in the vicinity of, the 
chiral restoration point, a mean field approach should make a better sense if 
one fluctuates around the "vector-manifestation (VM)" fixed point discovered 
by Harada and Yamawaki [10] of hidden local symmetry Lagrangian Wilsonian 
matched to QCD. Within the framework of hidden local symmetry (HLS) 
approach to hadron physics, there are two points at which the theory is reliably 
known. One is the zero-vacuum around which fluctuations can be computed 
with confidence if vector mesons are treated on the same footing as the pions. 
A suitable perturbation expansion treating the vector meson mass as a small 
parameter in the sense of the 1/N C expansion is found to give results in good 
agreement with experiments in matter- free space [10]. The other is the vector 
manifestation (VM) point at which in the chiral limit, the light-quark vector 
mesons become massless and the parameters of the Lagrangian controlling 
low-energy physics get fixed to a definite value and to which HLS flows as 
temperature, density or the number of flavors is dialled to the critical value 
for chiral restoration. If one is not too far from the VM point, it can be 
substantially advantageous to start from the VM fixed point at which the 
Lagrangian is known. We will call this "top-down" approach. An extreme case 
where this approach is successful even for a process taking place in matter- 
free space is the chiral doubler splitting of the D meson discussed in [11] 1 . 
Other examples are discussed in [12] where it has been suggested that certain 
processes in baryonic environment and/or particularly sensitive to the presence 
of vector degrees of freedom are much more effectively treated starting from 
the VM than from the zero-vacuum. This "top-down" approach is not yet 
fully formulated for the process in question. So we can only give a drastically 
simplified treatment below 



It is perhaps significant to note that in this process, even though the starting 
point is the VM, the tree contribution dominates, loop corrections making up only 
1/3 of the total. 



3 



In mean-field in the bottom-up approach on which we will be mainly focused, 
there are two main driving forces towards kaon condensation. One is the move- 
ment towards the restoration of the explicitly broken chiral symmetry T>kn i n 
the strangeness sector [13]. This is characterized in the calculations by the chi- 
ral symmetry parameter a^m s related to the KN sigma term. We shall update 
the calculations of strangeness condensation here of Thorsson et al. [1] who 
used three different values of a^m s . This parameter has now been calculated 
on the lattice [14], with central value a 3 m s = —231 MeV, and quoted accuracy 
of 3 — 4% which is only slightly greater in magnitude than the Thorsson et 
al. central value of —222 MeV which we shall use. Furthermore, Thorsson et 
al. used nuclear compression moduli of Kq = 180 and 240 MeV. We favor the 
value of 210 MeV for which = 1.5M . Kaon condensation sets in at 

n = 3.08n , where n is the density and n is nuclear matter density and the 
maximum neutron star mass is 1.5M . 

The second main driving force towards kaon condensation comes from the 
Weinberg- Tomozawa term treated in the mean field. This comes from vector 
mean fields between the nucleons and the K~ . Following [15], in mean field, 
the force mediated by the a;-meson exchange between a K~ and a nucleon can 
be related to that between two nucleons as 



V K -(cu) = - l -V N (1) 

where Vn is the vector mean field from nucleons in nuclear matter. The 1/3 
is easy to understand, because there is only one nonstrange antiquark in the 
K~ , whereas there are three nonstrange quarks in the nucleon. The p meson 
exchange gives 



V K -(p) = ±±V K -(u) (2) 

repulsive for neutrons, attractive for protons. In the Weinberg- Tomozawa 
term, the vector-meson propagator for symmetric nuclear matter is approxi- 
mated by zero momentum, 



Vir-(w) = -g^n, (3) 

where n is the nucleon density and x p (x n ) is the proton (neutron) fraction, 
the connection with the mean fields being by way of a KSRF-type relation 



m 2 v = 2Flg 2 v . (4) 
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In eqs.(3) and (4), F n is a parameter that appears in the HLS Lagrangian 
which is related to the pion decay constant f n 2 and gy is the vector (or 
flavor gauge) coupling. 

Two recent experimental developments provide crucial information that allows 
to determine the above mean field potentials in medium: 

(1) The "/„-" connected with the p- meson via the KSRF relation has been 
found in deeply bound pionic atoms to be substantially below its free 
space value [16] 3 , 

(i£(no)/F w ) 2 = 0.65 ±0.05 (5) 

extracted at tree order. Other determinations are reviewed by Brown and 
Rho [7]. 

(2) As what may be taken as the first direct and unambiguous verification of 
the scaling proposed in [17], the decrease in cu-meson mass with density 
has been measured as m*/m w ~ 0.84 at n = n [18]. This is somewhat 
larger than the 0.8, consistent with Eq. (5), which we shall use; we believe 
/* to be more accurately determined by the extensive data on deeply 
bound pionic atoms. 

As shown by Harada and Yamawaki [10], gy changes with scale, such that 
both g v and m v go to zero proportionally to the quark condensate (qq) near 
the VM fixed point as g v /m v constant. In fact there is some evidence that this 
behavior sets in from n ~ no upward [7,19,20]. However, there is no evidence 
that g v scales appreciably below n = Uq. Brown and Rho [7] have roughly 
summarized the situation by letting m v scale up to nuclear matter density 
as m v ( y no)/mv — 0.8, with gy constant. Beyond no gy/m v is taken to be 
constant. This is certainly a rough description, but it reproduces the main 
features of the Harada- Yamawaki scaling. 

With the scaling in F 2 only up to nuclear matter density no, we find the mean 
fields to be increased by a factor 



2 We make the distinction between and f^. It is which is an order parameter 
of chiral symmetry vanishing at the critical point in the chiral limit. Harada and 
Yamawaki show that f w = f 7T {Q 2 = 0) = F w +A where A is a quadratically divergent 
term that cancels F w at the critical point. Note that it is F n that we are concerned 
with here. 

3 We should note here that it is really F w that is determined in this experiment, not 
the order parameter f n . However at the mean field level and up to nuclear matter 
density, they are the same and hence the interpretation as "partial restoration of 
chiral symmetry" as measured in nuclei is correct. This aspect is frequently confused 
in the literature. 
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An astonishing new development has been the experimental discovery of deeply 
bound kaonic nuclear states [21]. The interpretation by Akaishi et al. [22] 
(see, for related discussions, Mares et al. [23] who studied the widths of 
deeply bound kaonic atoms and arrived at 100 ~ 200 MeV binding energy 
in these states) is of particular interest in connection with the Thorsson et 
al. [1] calculation and our present work is bringing it up to date. Although 
the deeply bound kaonic-nuclear states are small, essentially the strangeness 
(anti)nuggets suggested in many papers, the players in the nuggets, two neu- 
trons, one proton and a i^~-meson, are the same as those in the Thorsson et 
al. mean field approach, at the density beyond kaon condensation where one 
K~ is present for every two neutrons. We show that the Thorsson et al. mean 
field analysis gives substantial support to the Akaishi et al. analysis. In return, 
since medium dependent effects, essentially the dropping meson masses, have 
larger effects in the mean field than in the "nugget" calculations, experimental 
confirmation of the latter would give strong support to the work of Thorsson 
et al. and consequently the low M^| x of Brown and Bethe [2]. 

The plan of this note is as follows. In Sec. 2 we shall schematize the Thorsson et 
al. (bottom-up) calculation, carefully separating effects from scalar and vector 
mean fields in lowest order, showing how they get mixed up in higher order. 
We also discuss the range term, introduced in higher-order in [24,25] and show 
why the medium dependence in this effect has confounded investigators who 
obtain input parameters from i^~-nucleon scattering lengths. We then show 
in Sec. 3 how the mean field calculation is related to the strangeness "nugget" 
calculation of Akaishi et al. [22]. In Sec. 4, we discuss certain specific features 
associated with the VM point when the treatment is made top-down from the 
VM. In Sec. 5 we review some of the many other calculations of strangeness 
condensation and discuss why S~ hyperons do not seem to figure importantly 
in dense matter, or low — even if they do — they would not be expected to 
alter our M^g* appreciably. 

Although the detailed analysis of the evolution of compact binaries, double 
neutron stars and low-mass black-hole, neutron-star binaries, will be put off 
to a later paper [6], we point out the obvious consequences of a low M™g x 
mass in Sec. 6; namely that the giant progenitors of a binary neutron star 
must be very close (~ 4 %) in ZAMS mass. As more and more double binary 
neutron stars are observed, this can easily be checked. We shall point out 
that present observations already defy the statistics of the standard model of 
binary neutron star evolution, necessitating an alternative model. 

In Appendix, we construct a crystalized dense system comprised of tribaryon- 
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Fig. 1. Projection onto the a,K plane. The angular variable 6 represents fluctua- 
tions toward the kaon mean field. 

K~ bound states. We find that the K~ , s get unbound for densities exceeding 
2n . This indicates that in dense medium kaon condensation can set in starting 
from n ~ 2n . Since we started with the fixed proton fraction (1/3) in this 
approach, however, this critical density cannot be directly compared with the 
critical density for kaon condensation in the mean field approach. 



2 Strangeness Condensation 

In this section we sketch the principal arguments for our bottom-up approach. 



2.1 Scalar channel 

We first discuss the contribution towards the restoration of the explicitly bro- 
ken chiral symmetry, essentially considering V-spin [13] (see Fig. 1). This 
way of looking at kaon condensation does not capture all aspects of the phe- 
nomenon but it gives a clear conceptual picture of what is involved. 

The Hamiltonian for explicit chiral symmetry breaking is 



where the last expression is obtained for small fluctuation 9 in the direction 
of the K-axis. We then find 




(7) 
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Xkn(NN) \ 

fW K J 



(8) 



showing that the effective kaon mass m* K is strongly reduced by the explicit 
chiral symmetry breaking T, KN in the strangeness sector. 

As noted above, the explicit chiral symmetry breaking has been evaluated in 
lattice calculations by Dong et al. [14] giving T, KN = 362 MeV. Although in 
chiral perturbation theory the scalar nuclear density {4>ip) can be replaced by 
the vector density n = (ip^ip) since the above expression is in the highest order 
calculated, we prefer to keep the scalar density because at ~ 3no we may have 
gone over to constituent quark degrees of freedom. Looking at the behavior 
of the order parameter f*(n), we see that chiral symmetry restoration can 
only be declared 4 at n ~ 4n . Although one cannot take this result at its 
face value, we take this to suggest that the nucleon effective mass decreases 
linearly from n = to n c = Auq. This gives (ipip) / (ip^ip) for the proton as 
0.82, (ipifj) / (ip^if)) for neutron as 0.57 when we take 90 % neutrons and 10 % 
protons as Thorsson et al. [1] did at n — 3uq. So, the contribution of T<kn 
will be reduced by about 20% if we take the scalar density instead of vector 
density. In the arguments given below, instead of taking the scalar density 
self-consistently, we use vector density with smaller T^kn = 310 MeV, which 
correspond to a^m s = —222 MeV, for the comparison with the results of 
Thorsson et al.[l]. 

Now, one notable development following the Thorsson et al. work was the 
calculation of the range term [25] 



For self consistency the final in-medium co K - should be used here although 
this is at odds with the standard chiral counting that we are not adhering to. 
If we use the Thorsson et al. fi e = 219 MeV at n c = 3.08no for a^m s = —222 
MeV (remember that the K~ mass m* K - = uj* k - is equal to /i e at the phase 
transition.) then 



4 Whereas Nambu- Jona Lasinio theory gives chiral restoration at n c ~ 2uq [26] , the 
transition from nucleon to constituent quark degrees of freedom, which NJL starts 
with, should take another ~ 2no- With finite temperature the two transitions, one 
from nucleon to constituent quarks and then chiral restoration are discussed in [27] . 
It is possible to approach also n c ^ uq from the fixed point (at which m* — > in the 
chiral limit) since the role of Brown-Rho scaling is known (m*/m p ~ 1 — 0.2n/no) 
for n < no (see our later eq. (20)) and then the mass drops linearly in (qq)* for 
higher densities. 




(9) 
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Zf N = 0.93S XJV . (10) 
Without range term it follows from Eq. (8) that 

m* K = 333 MeV, (11) 
whereas with = 0.93£xjv 

m\ = 347 MeV. (12) 

Thus, with introduction of the range term there would be a 14 MeV correction 
upwards in the p e necessary for kaon condensation in Thorsson et al.[l]. 

The use of the full T, KN is highly preferable to the procedure used by most 
research workers who obtained from fitting the fT~-nucleon scattering 
lengths at threshold. They would have obtained 

^kn = 0.63£xjv (threshold scattering) (13) 
since ujk- = tuk- is assumed in this procedure. 5 

2.2 Vector channel 



We now turn to the Weinberg- Tomozawa term. As noted in Eq. (1), for uo- 
exchange the mean field isV k -{uj) = — (1/3) V/v, and for the p-meson V K - (p) = 
±(l/3)V K -(u). 

We calculate first with these mean fields. With the 90 % neutrons and 10 % 
protons at n c = 3.1n , the vector mean field contribution to cu K - would be 



V K -(u) = -\^\ + x P ) n ~ -126 MeV, ill; 



3 ml 2 V 2 

where g w = 3g p with g p ~ 5, and x n>p are the neutron and proton fractions. 
Putting this together with the m* K = 347 MeV of Eq. (12) gives us 



^K-( n — 3^o) = 221 MeV, (15) 

5 Indeed, for the pion the range term is ~ — LIS^jv and its introduction changes a 
quite appreciable attractive interaction into a slightly repulsive one. The fact that 
the range term in pion scattering is so well known makes it inexcusable to omit it 
in the i^T^-nucleon scattering. 
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essentially the same as Thorsson et al. [1]. However, from the medium depen- 
dence in hand, we increase the magnitude of the 126 MeV by 1.56 given by 
(6), giving an additional ~ 70 MeV drop in the K~ mass. As noted, 14 MeV 
of this is used up in the range term, so we are left with 56 MeV more binding 
than Thorsson et al. [1]. 

That the enhancement factor (6) is called for is clearly indicated in the struc- 
ture of nuclear matter. Were we to take the free pion decay constant f n ~ 93 
MeV in the K~ vector potential (3), it would give, at nuclear matter density, 



V K - = -V N ~ 58 MeV (16) 

which would mean that the cu-mean field for the nucleon would be V ~ 174 
MeV. This is much weaker than the Walecka mean fields which are employed 
at nuclear matter density n = n . They are more like ^(Walecka) ^ 270 MeV; 
i.e., ^ 50% higher. Thus, the phenomenology favors the mean field growing 
with density. 

In fact, the vector mean field is proportional to gv/m^ and a recent exper- 
imental study shows that m* decreases with density [18]. Much more study 
has been made of the density dependence of /* in deeply bound pionic atoms 
which arrive at 



(^p^-J = 0.65 ±0.05 (17) 

extracted at tree level [16]. 
A similar 



^-^0.8 (18) 
can be obtained [28] from the Gell-Mann-Oakes-Renner relation 6 



U m * = -Mqq) ( 19 ) 

under the assumption that m w and ifi do not scale with density, giving (/*/ f n ) 2 ■ 
(<?<?)*/ (<?<?); an d then using the model independent 



6 Note that here and in the next two equations, there is no difference between /„. 
and F v since we are going up to n ~ bq. 
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(qq) flml ' 

Thus by now we have ample evidence of scaling in m* and m*. 7 By taking 
/*/ f n ~ 0.8 in our formulae we assume them to drop from their bare value by 
20 % at nuclear matter density n . This gives our factor 1.56 increase in the 
vector mean fields. 



3 (Anti) Strange Nugget 



Recent experiments [21,29] and theoretical interpretation [22] involve the same 
players, e.g. two neutrons, one proton and a K~ in the strange tribaryon 
S' (3115) as in the neutron star matter. For example, in Table 4 of Thorsson 
et al. [1] which refers to a 3 m s = —222 MeV for n ~ 3.8no the proportion of 
K~ is x K - ~ 0.5, the K~ having replaced the electrons, so that the ratio of 
protons to neutrons is ~ 0.5 for overall charge neutrality. For this value of x 
the system is well into kaon condensation (with threshold u = n/n® = 3.08; 
note the similarity to mean density of S' (3115) of u = 3.1). Thus, if one cuts 
out a piece of the strangeness condensed star (which is on its way into a black 
hole) that contains two neutrons, it will have one proton and one K~ . 

First of all, in our description, the star should undergo strangeness condensa- 
tion at a lower value of u than the nugget. In the description of the star, mean 
fields are used; e.g., the Weinberg- Tomozawa term is 



9v 
3 m v 



(21) 



for the above composition of two neutrons for each proton. Note that since V ave 
is a mean field, it is evaluated for zero momentum transfer. This maximizes 
the medium effect in terms of the dropping vector meson mass my — > my. 

In fact, Brown et al. [30] have argued that in the GSI experiments with nucleon 
and kaon momenta \Pn\ ~ 444 MeV and \Pr\ ~ 322 MeV, form factors of 
fvip) ~ 0-82 must be employed in order to take into account the finite size 
of the nucleon. There is a partial decoupling of the vector interaction. Thus, 
the interactions in the nugget should be somewhat weaker than in the infinite 



The parametric masses of the to and p mesons should scale in the same way in 
medium in HLS theory if one assumes U(Nf) symmetry but the pole masses can 
differ a bit due to medium-dependent loop corrections. Here we are focusing on the 
former since we are working in the mean field approximation. 
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system, both because of the decoupling of the vector interaction and because 
medium effects are maximum for mean fields. 

The same interactions enter into the neutron star and strangeness nugget 
calculations, but the neutron star is well and truly strangeness-condensed, 
well beyond threshold for the composition at which the nugget is formed. 
This seems reasonable in terms of the stronger mean field interactions. 

Akaishi et al. [22] formulate their strangeness nugget problem in terms of po- 
tentials. Whereas the Weinberg- Tomozawa term, which is a vector potential, 
gives the main attraction between the two neutrons and the proton, we prefer 
to formulate the other main attraction as coming from the partial restoration 
of the explicitly broken chiral symmetry in the strangeness sector, essentially 
a scalar mean field. This would require something like a local density approx- 
imation in (ijjip) in the nugget problem, which would be of doubtful accuracy. 

We agree with Akaishi et al. [22] that the ~ 50 MeV difference between the 
binding energy (143 MeV) they obtained and the observed one (194 MeV) 
comes from the medium effects ( "Similar to the case of observed pionic bound 
states" [22] as we discussed in f v being replaced by /*) and it should be noted 
that in our description - which is different in spirit from theirs, we get about 
the same ~ 50 MeV as Akaishi et al. 

In other words, we believe that there is a mapping of the neutron-star problem 
onto the strangeness nugget problem, or vice versa, and that the discovery of 
the nugget strongly supports the former. 



4 a = 1 At and Away From the Fixed Point 

The behavior of the various quantities we work with should be better deter- 
mined at the fixed point of Harada and Yamawaki [10] where m* and g v both 
go to zero, linearly with (qq)* which goes to zero (as does /*) and the quan- 
tity a, which we now define — > 1. Mean field approximations are expected 
to get better the closer we are to the VM fixed point [10]. This suggests the 
top-down approach to kaon condensation as well as to the strangeness nugget 
problem. Here we present in the simplified form what we can say about kaon 
condensation. 

The low energy theorem of the hidden gauge theory valid to all orders of chiral 
perturbation theory [10] is 

< = ^K9v\ (22) 



12 



denning a in terms of the other variables we use. Earlier, when we talked about 
integrating out the p-mass, we used the KSRF-type relation 



m 2 p = 2Flg 2 v (23) 

which follows from vector dominance in free space, a being then equal to 2 
at the scale of m p . However, in turns out [31] that the vector dominance at 
a = 2 is on an unstable trajectory of RG flow of the HLS theory with no 
connection to the trajectory that leads to the Harada and Yamawaki vector 
manifestation and that the fact that in nature the vector dominance model 
seems to work in matter free space is merely an accident. In fact, apart from 
the pionic form factor at zero temperature and zero density, a near 1 provides 
a highly satisfactory phenomenology. In particular, a = 1 provides a good 
description for the coupling of the photon to the nucleon. [7] . Other evidences 
for a ~ 1 in nature including chiral doubling in D mesons are discussed in 
[12]. 8 Thus, we find it more appropriate to use 



a* 2 1 

% = — (24) 
with a* not far from 1, as compared with the 



9 A = — 2 (25) 

which pertained to matter-free space in which the vector meson was integrated 
out to get the Weinberg- Tomozawa relation. As can be seen, the vector mean 
field interaction is thus found to be increased by a factor of 2/a*, or a factor 
of 2 as the fixed point is reached. 

Our next objective is to estimate (in the chiral limit) at which density n c the 
fixed point is reached. We do this by finding out at which density m* goes to 
zero. 



Although initially m* decreases as y (qq)* following the scaling of F* (see 
argument above following Eq. (20)) the Harada and Yamawaki work shows 
that once it starts dropping it scales as (qq)*. (See the empirical verification 
of this in Koch and Brown [32] who showed that the entropy matched that in 
LGS if the meson masses were allowed to scale as (qq)*, which was referred to 



8 It is shown in [10] that at the matching scale A ~ 47r/ 7r ~ 1 GeV, HLS Lagrangian 
is given by a ~ 1.3. Since this Lagrangian is expected to hold in the large N c limit, 
one can think of a ~ 1.3 as representing the large -/V c limit. It turns out however 
that as far as phenomenology is concerned, a ~ 1 gives as good a fit as a ~ 1.3. 
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as "Nambu scaling.") As noted earlier, gy does not seem to scale up to nuclear 
matter density no, but then Nambu scaling sets in. Nambu scaling is times 
faster than the initial scaling of m* from n = to ~ no, which decreases m* 
by 20%. Thus, we believe in the interval n to 2no, m* will decrease \/2 times 
20%, or ~ 28%, and the same from 2no to 3no, and from 3no to nearly 4no 
where m* = in the chiral limit. Thus, the fixed point at n c is at n ^,4n . 

From our earlier argument that g v scales as m* for n > no, but up to n , gy 
remains constant, whereas m* scales, we find that Eq. (24) can be expressed 
as 



^--f— f (26) 

and we know that a* = 1 at n c , where it, together with m* and g v , has a fixed 
point. Compared with the matter-free expression Eq. (25), which is the KSRF 
relation, we see that 



[9y 2 / m f] fixed point _ _2_ _ g 13 ^ 

[gv 2 /m p 2 ] zero density 0.8 2 

Thus, the mean field felt by the K~ is increased by the factor 3.13 when 
the scaling of both F* and a are included, at the fixed point at n c , the final 
doubling coming from the scaling in a*. At the VM fixed point, the condensate 
is zero (in the chiral limit), so the scalar contribution from the rotation of Y>kn 
is gone. Thus seen from the VM point, the essential doubling of the attractive 
vector interaction replaces the attraction given by the rotating out of the 
Sxat term. In fact within the range of a relevant to the problem [12], say, 
1 ^ a 1.3, the coefficient of the density n needed to bring the m* K - down, 
is relatively constant up to n c . And furthermore the scenario is close to the 
bottom-up scenario discussed above. 

We emphasize that our above argument connects strangeness condensation 
with chiral restoration, in that with the large increase in vector mean fields in 
going to the fixed point, strangeness condensation certainly takes place at an 
n below n c . 

Following Eq. (15) it was noted that because of medium effects, we had an 
additional 56 MeV binding of the K~ as compared with Thorsson et al. Rather 
than a /i ejC = m* K -(n c ) ~ 219 MeV of these authors, we would find 221 — 56 
= 165 MeV. 

Now in our calculation about the fixed point, with the same 90% neutrons and 
10% protons at n c , we would have the factor of 3.13 times 126 MeV giving 394 
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MeV binding, just 100 MeV in magnitude less than the kaon mass in matter- 
free space. However, if we took a = 1.3 as suggested by large N c considerations, 
we would have 302 MeV binding, with jj, e c = tti^_ = 192 MeV. Thus, within 
the possible range of a that could be arrived at by the RG flow from the fixed 
point (we believe ~ 1.3 to be maximal), we find the \i efi of Thorsson et al., 
corrected by medium effect, to be midway in the range of // e , c 's that could be 
reached starting from the fixed point. A similar attraction was obtained for 
rip/n ~ 0.25 by Tsushima et al. [33] in a different approach that takes into 
account the scaling behavior. 

We see from Thorsson et al., for a = 1 — 1.3, the necessary electron chemical 
potential is easily reachable by n = 3no so strangeness condensation seems 
assured by 3n . In fact, the mass of the Hulse- Taylor pulsar is 1.44M and 
since it must be stabilized, it will give a lower bound on n c . 

The Kaplan- Nelson term involving Y^kn that is important in the bottom- up 
approach to kaon condensation in an infinite matter would be difficult to 
include reliably in the strangeness nugget calculation, so it would be far more 
advantageous to carry it out top-down from the VM point. Since they do not 
deal with mean fields, but rather with relatively high vector meson momenta 
in 



9v 2 

mtf + q* (28) 

they will not benefit as much as the kaon condensation calculation from the 
change in a from 2 in free space to ~ 1 in medium. However, the factor 1.56 for 
(F v /F*) 2 remains. Thus, we might expect the binding there to be ~ 1.56 x 126 
MeV =147 MeV. In fact 126 MeV was the original vector mean field for 90% 
neutrons and 10% protons whereas Akaishi et al. have two neutrons for each 
proton, so the number would be 160 MeV at n — 3fiQ. Hence the binding 
would be 1.56 x 160 MeV = 250 MeV. 



5 The Lack of E" Role 



The mass of the S~pn system lies at ~ 3075 MeV, well below that of the 
S' (3115), and one might think that this system plays a role in strangelets. 
Similarly, in neutron stars, the S~ will replace both a nucleon and an electron. 
This consideration led many investigators to consider £~ condensation, most 
recently by Kolomeitsev and Voskresensky [34]. 

We do not think that the £~ obstructs our principal argument. Batty et al.[35] 
noted that "for K~ atoms, the fitted potential becomes repulsive inside the 
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nucleons, implying that E hyperons generally do not bind in nuclei." One can 
see from Fig. 4 in Brown et al. [36] that in the density dependencies favored 
by Batty et al. [35], the S~ optical potentials are highly repulsive. Further- 
more there are now direct experimental nuclear data (to be distinguished from 
atomic data) that show that the repulsion increases with the neutron excess 
[37,38]. A DWIA analysis [39] of these (7r~, K + ) data confirms the repulsion 
which was first predicted from S~ atoms by Batty et al. [40] and extended by 
Mares et al. [41] in a relativistic mean field approach. 

In Kolomeitsev and Voskresensky [34], although the S-wave K~ condensate 
occurs only at p c ^ 4p in neutron star matter, it is preceded by a long, 
essentially mixed phase of hyperon condensation which begins at u ~ 2. During 
a mixed phase the pressure is constant; even with Gibbs construction the 
gradient in pressure is low, so that the force, given by gradient in pressure, is 
low. Thus, the neutron star should compress substantially in this phase and 
it seems unlikely that the maximum mass will be greater than our 1.5M Q . 



6 Astrophysical Consequences 

In this Section we briefly touch on astrophysical implications relegating details 
to a future publication [6]. The important astrophysical consequence of the 
low-mass MjSfg x is that the standard scenario for binary neutron star evolution 
[42] does not result in a double neutron star, but in a black-hole, neutron-star 
binary. In this scenario, after the first born neutron star is formed, it goes into 
common envelope evolution with the companion giant as the latter expands in 
red giant stage. During this common envelope evolution it accretes a substan- 
tial amount of matter from the hydrogen envelope of the giant companion, as 
it removes this envelope. Bethe and Brown [43] estimated this amount to be 
~ 1M for a 1.4M neutron star, whereas with more accurate calculation [44] 
found ~ 3/4M for this mass neutron star. The accretion is ~ O.5M for a 
1.1 — 1.2M neutron star. Obviously these will be sent into black holes and 
the result will be a black-hole, neutron star binary. 

As shown in [45], there is a special way in which the neutron star evolution 
in common envelope could be avoided. If the two giant progenitors are 4% 
different in mass, they will expand in red giant and burn helium at the same 
time. Their hydrogen envelopes are then removed in the double helium star 
common envelope evolution. There is no time for the hydrogen to cross the 
helium molecular weight barrier [46] so the helium star remains very close in 
mass. The result is that the two resulting neutron stars have very nearly equal 
masses. 

In the Table 1 we show the presently measured double neutron star binary 
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Table 1 

Compilation of the compact objects in binaries by Lattimer and Prakash (2004). 
References are given in their paper. We have added, following the comma, the recent 
measurement of Van der Meer et al.[50]. 



Object 


Mass (M ) 


Object 


Mass (M Q ) 


X-ray Binaries 








4U1700-37 


—0.27 


Vela X-l 


1.86+n{6 

—0.1b 


Cyg X-l 


7 o+0.23 
—0.26 


4U1538-52 


0.96+°^ 

— U.lO 


SMC X-l 


1. 17 + n]L 1.05+0.09 
— u. id ' 


XTE J2123-058 


1 ro+0.30 
X . (JO n /I o 
— U.4z 


LMC X-4 


1.47+S-^, 1.31+0.14 


Her X-l 


— U.lo 


Cen X-3 


1.09+H9, 1.24+0.24 


2A 1822-371 


> 0.73 


Neutron Star - 


Neutron Star Binaries 






1518+49 


—0.44 


1518+49 companion 


""—0.11 


1534+12 


1 Qqqo+0.0010 
—0.0010 


1534+12 companion 


1.3452+°-™}° 

—0.0010 


1913+16 


1.4408+Hons 

—U.UUUo 


1913+16 companion 


1 0070+0.0003 

— U.UUUo 


2127+11C 


1.349+g-gS 

— U.U4U 


2127+11C companion 


1 

— U.U4U 


J0737-3039A 


i OO7+0.005 
— U.UUo 


J0737-3039B 


1-250+S-SSi 

—U.UUo 


J1756-2251 


—0.03 


J1756— 2251 companion 


1 1Q+0.03 

1 . lO n no 

— U.U2 


Neutron Star - 


White Dwarf Binaries 






B2303+46 


1 OO + 0.06 

J-^-o.io 


J1012+5307 


1 68+^2 

uu -0.22 


ti 71 q_i_0747 


i tr 4 +0.007 
-L-^-O.OOS 


B1802-07 


1 26+ - 08 


B1855+09 


i c 7 +0.12 


J0621+1002 


x - ' u -0.29 


J0751+1807 


9 9n +0.20 
z - zu -0.20 


J0437-4715 


1 co+0.18 
l-O8_ 18 


J1141-6545 


1 30+ 002 


J1045-4509 


< 1.48 


J1804-2718 


< 1.70 


J2019+2425 


< 1.51 


Neutron Star - 


Main Sequence Binaries 




J0045-7319 


i - O5 -0.34 







masses. In the Hulse- Taylor pulsar 1913+16 the pulsar mass is 1.44M , the 
most massive neutron star in the double neutron star binaries. It has been 
evolved from an ~ 2OM giant [47] . Since the number of giants go as 



N(M) oc ^ (29) 
according to the Salpeter mass function, if there were no correlation in pulsar 
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and companion masses, its companion would be most likely to come from an 
M ~ 1OM giant, with neutron star mass ~ 1.2M Q . 

As can be seen from the accurately measured 1534+12, 1913+16, and 2127+1 1C 
there is a close correlation in the masses of the pulsar and companion. 

In the case of the lower-mass binaries, the double pulsar J0737— 3039A,B and 
J1756— 221 equality of masses has been disturbed by the helium star evolu- 
tion. In the chain of evolutions, the first born neutron star has a helium star 
companion, and low-mass helium stars expand in their own red giant stage. 
In so doing they transfer ~ 0.1 — O.2M to the neutron star, but since these 
are of low mass, this is not enough to send them into a black hole. 

Since substantial mass is transferred in the neutron-star, hydrogen-envelope 
common envelope evolution for the more massive stars, the above correlation 
in neutron star masses within a binary does not give an accurate limit on 
m^g x . However, Brown's [45] special double helium star scenario is an order 
of magnitude less likely than the standard scenario, giving ~ 10 times fewer 
double neutron star binaries. This translates into 10 times more black-hole, 
neutron-star binaries than double neutron-star binaries, which, with the higher 
masses of the former because of accretion, should give a factor of ~ 20 increase 
in the gravitational waves from merging to be detected by LIGO [43]. 

Since Table 1 contains three masses, those of 4U 1700—37, Vela X-l and 
J0751 + 1807 which exceed our 1.5M maximum neutron star mass, we should 
comment briefly. 

4U 1700—37: Although this compact object has the same accretion history 
as the other high-mass X-ray binaries, it doesn't pulse like the others. Brown, 
Weingartner and Wijers [5] evolve the compact object as a low-mass black 
hole. 

Vela X-l : J. van Paradijs et al. [48] pointed out that in this binary with 
floppy B-star companion, the apparent velocity can in some cases increase 
by up to 30% (from the surface elements of the companion swinging around 
faster than the center of mass) "thereby increasing the apparent mass of the 
compact object by approximately the same amount". In any case, Barziv et 
al.[49] from which the Vela X-l neutron star mass in our table comes, say 
"The best value of the mass of Vela X-l is 1.86M Q . Unfortunately, no firm 
constraints on the equation of state are possible, since systematic deviations 
in the radial-velocity curve do not allow us to exclude a mass around 1.4M Q 
as found for the other neutron stars." 

J0751+1807: We are unable to give a critical discussion of this binary because 
thus far only preliminary results have been published. The sum of neutron-star 
and white-dwarf mass are determined by the period change from gravitational 
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radiation. There is an indication of a Shapiro shift and a white dwarf mass of 
0.188 ± 0.012M Q is arrived at. Details of this are not given (of course it could 
be obtained from the measurement of the Shapiro shift, but none is given.) 
With a white dwarf mass of O.24M , the neutron star mass would come down 
to 1.5M . 

We should mention that the existence of a 2.2M neutron star would wreak 
havoc with our scenario of nearly equal masses in the more massive of the 
binary pulsars. The most abundant first-born neutron stars of mass 1.1 — 
1.2M Q from ~ ZAMS 1OM giant progenitors would accrete ~ O.5M in 
hydrogen red giant phase and another 0.1 — O.2M in helium red giant. They 
would then end up with pulsar mass ~ 1.8M , with most likely low-mass 
1.1 — 1.2M companion. Such massive pulsars should be strong in the radio 
and should predominate in number. None such have been observed. 



7 Discussion 

Although our discussion based on the bottom-up approach supports the Akaishi 
et al. scenario of strangeness nuggets, we have emphasized the Kaplan and 
Nelson [51] scalar attraction from the movement towards restoration of the 
explicitly broken chiral symmetry in the strange sector, Eqs. (7) and (8). Of 
course, this is more straightforward for us to handle in our mean field ap- 
proach, than in the strangeness nugget, although the ratio of (ipip) to vector 
density n is model dependent. On the other hand, the vector mean fields such 
as Vk-{oj) = — \Vn of Eq. (1) are straightforward, the K~ having one non- 
strange antiquark, and the constituent quark model should be adequate at 
n = 3no to describe the | ratio of Vk- to Vn- The vector mean field Vn at 
n = n has been much used in Walecka theory in nuclear physics. A value of 
V/v(no) = 270 MeV is quite modest in nuclear physics. It should be noted that 
this value of 3Vk-(ou) only reaches this after the factor (F^/F*) 2 is included. 
There is a clear indication, not least of all directly from experiment [18] that a 
sliding vacuum; i.e., Brown- Rho scaling must be operative. The vector meson 
masses continue to drop, until they go to zero in the chiral limit, but as the 
scale (density in this case) increases, the coupling g v also decreases. Near the 
fixed point {gy/rriy) 2 goes as a constant. We have handled the scaling crudely 
by assuming that gy does not scale up to n , for which there is substantial 
basis [7], and then neglect any further scaling in g v /m v above n . This gives 
the factor (F w /F*) 2 = 1.56 in the medium dependence we used. In the above 
way we achieve a rough consistency with the Harada and Yamawaki [10] hid- 
den local symmetry. Since all of our important results here refer to n > no, 
for these the sliding in coupling cancels that in meson mass, and we can carry 
out the calculations with an effective /* such that (F n /F*) 2 = 1.56 replacing 
This was also foreseen in the Akaishi et al. work. In short, Brown-Rho 
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scaling essentially increases the zero-density parameters of mean fields up to 
the couplings which have already been used for a long time in Walecka type 
mean field approaches. 

We suggested that mean field would be much more reliable when treated start- 
ing from the VM fixed point and sketched a highly simplified calculation that 
reproduces the bottom-up mean field calculation for kaon condensation. When 
further developed, the top-down approach will prove to be more powerful for 
the strangeness nugget problem. 

It must be admitted that standard chiral perturbation theory calculations 
fluctuating from the zero- vacuum have not been carried out to the order nec- 
essary to get from the input lowest-order f n taken to be equal to fx in the 
necessary places. It is not clear to us that one can obtain a reliable result from 
such an endeavor. It seems that such a calculation must begin in-medium, 
from F* ~ 0.8F„-. Thus, the in-medium would not be expected to be larger 
than the f n we have used in m* K /rriK- 

An interesting calculation in chiral perturbation theory was carried out by 
Waas and Weise [52] and did go to high enough order to include the range term 
(they did not however employ the medium modified /*). They used T,kn = 230 
MeV, substantially smaller than that found in LGS. Waas and Weise found 
that iterating the Weinberg- Tomozawa term increased the attraction, lowering 
m^/rriK appreciably. In this way, their approach was midway between our 
mean field and the Akaishi et al. approach, which iterated the attraction by 
solving in a potential. 

In our discussion of the melting of the soft glue in LGS [53] we found that 
nucleons changed over to constituent quarks at T ~ 120 MeV whereas the 
constituent quarks became current quarks, having lost their dynamically gen- 
erated masses, at T c (unquenched) = 175 MeV. Although we do not know the 
value of n c , and as we noted earlier, in Nambu Jona-Lasinio (ipip) doesn't re- 
ally go to zero until well above the n c for zero bare quark masses; we suggest 
that at u = n/riQ ~ 3 constituent quarks may be more appropriate variables 
than nucleons, although they will still be held together by the soft glue in 
nucleons to some extent. 

Thus, although fx may enter in the Kaplan-Nelson term Eq. (8), the vector 
mean fields are connected only to the constituent u quark in the K~, and 
we believe that all of these calculations can be handled at the constituent 
quark level, but with Brown-Rho scaling of the constituent quark masses. 
This suggests working around the VM fixed point as stressed in [12]. Whereas 
the lowest-order Weinberg- Tomozawa term brings m* K /mK down to 0.65 in 
nuclear matter (it would not be brought down so far in neutron stars) with 
a factor of 1.56 it would be brought down more than half way to zero. The 
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second order terms are necessarily attractive, and will be increased by a factor 
of (F V /F*) A ~ 2.4. There may be saturation in the complete coupled chan- 
nels calculation, which includes Pauli blocking, but the contribution of the 
Weinberg- Tomozawa terms will obviously be substantially increased by the 
medium dependence. On the other hand, the (ipip) for constituent quarks may 
be smaller than the value Thorsson et al. [1] calculated for nucleons. We be- 
lieve that the additional attraction from medium effects on the vector mean 
fields would be adequate to compensate for this. 

The strangeness nugget presents a remarkable scenario in which once the K~ 
meson is introduced into the three remaining nucleons in 4 He, after a pro- 
ton has been kicked out, the three nucleons collapse around the K~ into the 
small nugget system. What else can they do ? All interactions, scalar from 
the Kaplan-Nelson term and vector mean fields from Weinberg- Tomozawa 
terms are highly attractive, the latter having only a small isospin dependence 
through the p mean field, so the system collapses down to an average den- 
sity n ~ 3riQ. (It is however difficult to imagine how this can be achieved in 
standard chiral perturbative approaches.) The K~ is in a potential about 600 
MeV deep at the center, with binding energy of ~ 250 MeV, even more after 
strangeness condensation in our neutron stars. Luckily the collapse is halted 
in the nugget formation, chiefly by the finite sizes of the interacting objects. 
In the case of neutron stars, there is collapse into black holes and we will not 
"see" the black hole, neutron-star binaries until gravitational waves from their 
mergers are detected. 

One issue we have not addressed in this paper is whether the collapse to black 
holes occurs without the kaon condensed state going into a color supercon- 
cuting quark matter or after such a transition. At the moment one can say 
practically nothing about this matter. Should the collapse occur from a quark 
matter rather than from a kaon-condensed state, then the question as to how 
the kaon-condensed state which we predict must occur before the chiral tran- 
sition and which may not be a normal Fermi liquid state develops the pairing 
instability that turns the matter to a (color) superconducting state discussed 
in the literature. As far as we know this question has not been studied in the 
literature. 



8 Conclusions 

We have argued that strangeness condensation in neutron stars can be related 
to the strangeness nuggets found experimentally by Suzuki et al. [21,29]. In 
both cases, binding of the K~ meson is responsible for the phenomenon. In 
the neutron star matter the applicable mean field method gives ~ 50 MeV 
(or ~ 20%) more binding to the K~ than in the nugget, where m* K - ~ \tuk- 
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In both cases, the medium dependence in the vector meson mass m* and m* 
provide ~ 50 MeV of the binding. Our arguments are admittedly far from 
rigorous. If this relation can be put on a more rigorous basis, it will supply 
a firm support for out limit of 1.5M as maximum neutron star mass. It is 
remarkable that in Table 1 out of ~ 40 measured masses of compact objects 
at most two or three of the masses seem to violate our upper limit. These and 
other astrophysical issues will be discussed in depth in a paper in preparation 
[6]- 



Dedication 

We would like to dedicate this paper to Hans Bethe who fearlessly coauthored 
"A Scenario for a large number of low-mass black holes in the Galaxy" [2], the 
large number resulting from the low maximum neutron star mass of 1.5M . 
Gerry Brown was explaining the idea of kaon condensation to Hans Bethe 
on a Saturday morning walk on a trail near Santa Barbara, California. Hans 
immediately understood it, "You mean that you squeeze electrons into Ja- 
mesons ?" 
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A Appendix: Dense system with tribaryon-A' crystals 

In this Appendix, we construct a crystalized dense system comprised of tribaryon- 
K~ bound states. The assumptions we make in the construction of the crystal 
potential are: 

(1) The effective potential of the single tribaryon-A" (ppn-A', Isospin=l) pro- 
posed by Akaishi et al. [22] is used; 

V{r) = (V + W )exp[-(r/a) 2 ] (A.l) 

with V = -702 MeV, W = -13 MeV and a=0.923 fm. By solving the 
Klein-Gordon equation, 
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nj no 


1 DQ 1 ^7 1 7i 1 Q7 9 (19 

i.uy i.oi ±.y/ z.uz 


B.E.(MeV) 
rms Radius(fm) 


-217.7 -238.5 -275.4 -304.3 -310.7 
-136.6 -188.3 -258.1 -302.2 -311.9 
0.816 0.842 0.862 0.868 unbound 







Table A.l 

Binding Energy and rms radius of a single tribaryon- K system in medium for various 
densities. Beyond 2no, the bound tribaryon- if - system "dissolves" indicating that 
the the crystalized system is unstable. 

((E-V(r)) 2 -m 2 K -p 2 }^ = 0, (A.2) 

Akaishi et al. obtained the binging energy and rms radius of the tribaryon- 
K as 194 MeV and 0.74 fm, respectively. In this estimate, the enhanced 
core energy of the tribaryon from the core shrinkage, 5E ~ 110 MeV, 
is assumed. Akaishi at el. show the average density of the tribaryon-.^" 
bound system to be n ~ 3.1no- 

(2) We put these tribaryon- K~ systems in a cubic lattice by assuming that 
the potential generated by the tribaryon is the same as that given by 
Akaishi et al. 

(3) In a dense medium, a tribaryon- .K" bound system is affected by the sur- 
rounding tribaryon cores. In order to take this effect into account, we 
make an effective potential of a single site by making the average over all 
solid angle. Finally we have an effective potential with spherical symme- 
try. 

With the above assumptions, we obtain the binding energies and the rms radii 
of the tribaryon- K~ system. In order to take into account the density effects, 
we change the distances between these tribaryon-fT~ systems. 

In Figure A.l, we see that the minimum and the threshold of the effective 
potential of a tribaryon- K~ system are lowered as the baryon density increases 
due to the overlap of the potential given by the neighboring systems. As one 
can see from the figure, when the binding energy is equal to the maximum 
value of the potential, the bound state of the tribaryon-^ system dissolves. 

As in Table A, the bound state of the tribaryon- .K" disappears at n/no — 
2.0. This indicates that for densities greater than 2n , the system becomes 
continuous, losing the crystal structure. We interpret this as kaon condensation 
taking place from 2n . Since we started with the fixed proton fraction (1/3) 
in this approach, the density at which the crystal dissolves cannot be directly 
compared with the critical density for kaon condensation in the mean field 
approach. However, our results provide a way to link tribaryon-i^~ bound 
states (nuggets) to kaon condensation. 
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Fig. A.l. Effective potentials for various baryon densities. 
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